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ABSTRACT 


The  dynamic  analysis  of  "the  wall  of  a fluid-filled 
unstiffened  nuclear  containment  vessel,  to  the  fluid  pres 
sure  exerted  on  it  when  the  relief  valve  discharge  piping 
is  cleared,  is  extended  into  the  plastic  range  using  two 
versions  of  an  elastic-plastic  shell  theory. 


M.>  ... 


1.  INTRODUCTION 

In  an  earlier  study,  (Ref.  [l]),  a method  was  presented 
for  obtaining  the  dynamic  elastic  response  of  a fluid-filled 
nuclear  containment  vessel  when  the  relief  valve  discharge 
piping  is  cleared.  In  this  thesis,  the  method  is  extended 
into  the  plastic  range,  using  modified  versions  of  a plas- 
ticity theory  for  shells  introduced  by  Bieniek  and  Fumaro 
(Refs.  [2], [3]).  Alternatively,  this  investigation  could 
have  been  conducted  without  using  a shell  theory.  Instead, 
the  field  equations  of  classical  elastic-plastic  theory 
could  be  applied  at  shell  layers  and  then  integrated  through 
the  thickness  (see,  e.g..  Refs.  [4]  to  [8]).  This  latter 
procedure,  however,  may  require  prohibitively  large  comput- 
er storage  capacity. 

The  water-filled  circular  cylindrical  shell,  shown  in 
Fig.  1,  is  the  model  used  for  the  lower  portion  of  a nucle- 
ar containment  vessel.  It  is  acted  upon  by  a cylindrical 
( axisymmetric  ) wave  whose  time  history  is  shown  in  Fig.  2 
( with  magnitude  normalized  to  1 psi)  and  whose  variation 
with  depth  is  illustrated  in  Fig.  3.  This  approximates 
(see  Ref.  [9] ) the  pressure  exerted  on  the  vessel  wall  when 
the  relief  valve  is  cleared.  The  spatial  and  temporal  var- 
iation of  displacement,  moment,  and  hoop  stress  in  the  cy- 
lindrical steel  wall  is  to  be  determined. 


2.  FORMULATION  OF  THE  PROBLEM 


2.1  Fluid  Equations 

Assuming  linear  compressibility,  the  pressure  p(r,z,t) 
in  the  contained  water,  and  its  (vector)  velocity  u(r,z,t) 
satisfy  the  equations  (Ref.  [lO]) 


3a  _ 

e-jt — v»  <2> 

in  which tp  is  the  density  of  water j c is  the  velocity 
of  sound  in  water » and  V is  the  (vector)  gradient  oper- 
ator, and  the  dot  denotes  scalar  multiplication. 

The  continuity  requirement,  that  the  radial  fluid  velo- 
city at  the  shell  surface  be  equal  to  that  of  the  shell, 
may  be  written  as 


-ff-  - u(rt)  (3) 

in  which  u(r,z,t)  is  the  radial  component  of  the  fluid  ve- 
locity u,  r is  the  radial  coordinate  as  in  Fig.  1,  and 
w is  the  radial  displacement  of  the  shell,  assumed  positive 
outward . 
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The  boundary  conditions  to  be  satisfied  arei 


w(l,t)  » 0 


^(l.t)  - 0. 

oz 


(4) 


(5) 


Equations  (4)  and  (5)  assume  that  the  shell  bottom  is  fix- 
ed. 1,  as  shown  in  Pig.  1 , is  the  depth  of  the  shell 
wall. 
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(0,t)  = 0 


(6) 


(7) 


Equations  (6)  and  (7)  assume  the  shell  top  to  be  free. 


u(r1#z#t)  * 0 (0  ^ z ^L) 


(8) 


Equation  (8)  assumes  the  pedestal  (rigid  cylindrical  sur- 
face of  Fig.  1)  to  be  rigid,  so  that  the  radial  fluid  vel- 
ocity is  zero  there.  L is  the  depth  of  the  flat  rigid 
bottom  from  the  free  water  surface,  as  shown  in  Fig.  1. 


v(r,L,t ) =0  (r^r<r2). 


(9) 


Equation  (9)  states  that  the  flat  portion  of  the  vessel's 
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bottom  (see  Fig.  1)  is  rigid,  requiring  that  the  axial  com- 
ponent v(r,z,t)  of  the  fluid  velocity  u vanish  there. 

[u  si„«  + V cob ,=1+^r^-r)taiiK * 0 (10) 

Equation  (10)  states  that  the  sloping  bottom  of  the  con- 
tainment vessel  (see  Fig.  1)  is  rigid,  so  that  the  compo- 
nent of  fluid  velocity  normal  to  it  is  zero. 

p(r.O.t)  « 0 (r^  r ^ r3) . (11) 

Equation  (11)  states  that  the  fluid  pressure  on  the  free 
surface  is  zero. 

Initial  rest  conditions  require  that 


w(r,0) 

= 0 

(12) 

dw 

— (z,0) 
dt 

= 0 

(13) 

p(r, z,0) 

= 0 

(14) 

u(r,z,0) 

= 0 

(15) 

Note  that  Eqns  (2),  (3),  (11),  (12)  and  (14)  constrain 
the  shell  displacement  at  the  fluid  surface  to  be  zero,  i.e. 


w(0,t)  * 0 


(16) 
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2. 2 Shell  Equations 

From  the  radial  component  of  the  shell  displacement, 
w(z,t),  assumed  positive  ou+^ard  as  shown  in  Fig.  1,  the 
circumferal  strain  € and  longitudinal  curvature  ^ of  the 
shell  are  obtained  as 


w 

r- 


(17) 


a7 


(18) 


The  corresponding  stress  resultants  are,  respectively, 
the  circumferential  force  per  unit  of  length,  N,  and  lon- 
gitudinal moment  per  unit  of  length,  M,  which  satisfy  the 
stress  equation  of  motion  (see,  e.g..  Ref.  M ) 


d2M 


N 


32 


■—?  ♦ - PT(z,t)  + p(rvt)  - p h 

dz*  r3  J s 


(19) 


in  which  E is  Young's  modulus,  Pj  is  the  incident  pres- 
sure of  Figs.  2 and  3 caused  by  the  relief  valve  clearing, 
pg  is  the  mass  density  of  the  shell,  and  h is  the  shell 
thickness. 

The  total  strains  of  Eqs.  (17)  and  (18)  will  be  ob- 
tained as  a superposition  of  elastic  and  plastic  components « 
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Using  matrix  notation, 


Eqs.  (20)  and  (21)  become 


e = e^  + e** 


(23) 


2.2.1  Elastic  Range 


The  elastic  stress  strain  relations  are  (Ref.  [ll]  ) 


fV 

12(l-y^) 


N = Eht 


(25) 


in  which  y is  Poisson’s  ratio,  E is  Young’s  modulus,  and 
h is  the  shell  thickness.  Letting 


«•(:] 


12(l-y  ; 


(27) 


be  the  stress  and  elastic  modulus  matrix,  respectively,  Eqs, 
(24)  and  (25)  may  be  written  in  matrix  form  as 


s = E eB  * E (e  - eP) 


Substituting  Eqs.  (17),  (18),  (24)  and  (25)  into  the 
stress  equation  of  motion,  the  elastic  displacement  equa- 
tion of  motion  is  obtained  as 


__  £5  ♦ - - . Pl  ♦ P<r3>  - p,h  —2 


2.2.2  Elasto-Plastic  Range 


The  shell  will  be  assumed  to  yield  when  the  yield  con- 


dition 


F0  = IN  + XM  + 2 rMN  = 1 > 


in  which 
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(33) 


2M0N0 


is  satisfied.  This  relation  is  obtained  from  that  general 
three  dimensional  (condition)  used  by  Bieniek  and  Funaro 
(Ref.  [2]  ),  for  the  axi symmetric  condition  of  the  problem 
being  considered,  when  Poisson's  ratio  effects  (i.e.,  in- 
duced circumferential  moments)  are  neglected.  In  Eqs. 

(3D  - (33). 


No  = *oh 


Mo  = 


in  which  0^  is  the  uniaxial  yield  stress*  NQ  and  MQ  are, 
respectively,  the  axial  force  at  which  yielding  occurs  in 
the  absence  of  bending,  and  bending  moment  at  which  yield- 
ing occurs  in  the  absence  of  axial  force,  in  a beam. 

If  circumferential  moments  are  again  ignored,  the  li- 
mit function  suggested  in  Ref.  [2]  is 


PL  " * 


IM  ♦ 


M 3fr 


(36) 
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In  Pig.  4,  the  yield  and  limit  surfaces  corresponding  to 
Eqs.  (30)  and  (36)  are  displayed. 

In  conjunction  with  the  yield  and  limit  functions  of 
Eqs.  (30)  and  (36),  it  is  proposed  to  use  a loading  function 


F 


♦4* 


2 

3 JT 


(37) 


in  which 


# M2  - 2MM* 


In  Reference  [2]  , 


(M  - M*) 


2 


(36) 


(39) 


is  used.  For  stress  points  near  the  axis  M = 0 the  result- 
ing loading  function  can  lead  to  a violation  of  Drucker's 
second  postulate  (Ref.  [12]  ) 


(40)* 


Drucker's  first  postulate  (Ref.  [12]  ) 


(s  " Sq)  eP  ^ 0, 


(41) 


* In  the  Appendix,  an  alternate  procedure  for  overcoming 
this  difficulty  is  considered. 


in  which  sQ  is  any  stress  point  within  the  loading  function 
F=l,  is  satisfied  using  either  expression  for  I*  since  the 
resulting  loading  surfaces  are  convex  and  an  associated 
flow  rule  will  be  used. 

As  in  Ref.  [2]  , it  is  proposed  that  the  hardening 

* . 

parameter  M increase  during  plastic  loading  according  to 
the  rule 


in  which  MQ  and  kQ  are  related  by  Eq.  (24), 


3(N/Nq) 


is  the  absolute  value  of  the  gradient  of  the  loading  sur- 
face in  nondimens ional  space,  and 


d(M/ta0) 


is  that  part  of  G contributed  by  bending. 

The  resulting  loading  function,  shown  in  Fig.  4,  ex- 
hibits features  of  both  kinematic  and  isotropic  hardening. 
Wherever  F=1  lies  outside  FL=1,  FL=1  is  used  to  determine 
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Plastic  loading  is  defined  to  occur  when  Eq,  (37)  is 
satisfied  and 


dP  . 3F  . 

SF  = — N + — M>0  (45) 

8N  9m 


in  which  dots  denote  differentiation  with  respect  to  time. 
(Unloading,  which  is  elastic,  is  characterized  by  Sf<0, 
while  for  neutral  loading,  during  which  the  stress  point 
moves  on  the  loading  surface  but  no  hardening  occurs,  §F=0). 
Letting 


3F 

3N 


9s 


BF 

w 


(46) 


The  loading  condition  of  Eq.  (45)  may  be  written  as 


Sf 


3f 

9s 


s > 0 


(4?) 


in  which  the  symbol  denotes  transpose. 

The  plastic  strain  increment  is  obtained  from  the  load- 
ing function  by  an  associated  flow  rule 


= X 


9F 


(48) 


9s 
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is  the  elastic-plastic  tangent  stiffness. 

It  should  be  noted  that,  substituting  Eq.  (52)  into 
Eq.  (40),  Drucker's  second  postulate  becomes  equivalent  to 
the  requirement  that  D be  positive  definite. 
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3.  METHOD  OF  SOLUTION 

The  finite  element  method  utilized  for  the  fluid  equa- 
tion and  the  elastic  phase  of  the  shell  motion  are  identi- 
cal to  those  described  in  Ref.  [l]  and  will  not  be  repeated 
here.  The  two-dimensional  grid  and  molecule  used  for  the 
fluid  equation,  like  the  one  displayed  in  Fig.  5,  and  the 
one  dimensional  grid  used  for  the  elastic  phase  of  the  shell 
motion,  shown  in  Fig.  6,  are  also  used  for  the  plastic  phase. 

In  what  follows  the  iterative  computational  scheme 
utilized  in  the  (nonlinear)  plastic  phase  is  outlined.  In 
all  cases  where  equations  are  referenced,  the  finite  diffe- 
rence form  of  the  equation  was  used. 

Assume  that  at  a time  t^,  when  the  incident  pressure 
Pl(ti)  is  acting,  the  stress  s(ti)  strain  £(t^)  and 
displacement  w(t^)  in  the  shell,  and  the  velocity  u(t^) 
and  induced  pressure  p(t^)  corresponding  to  a grid  point 
are  known.  In  this  notation  and  that  which  follows  in  this 
section,  the  spatial  dependence,  usually  expressed  with  sub- 
scripts in  finite  difference  form,  is  supressed.  Let  M*(ti) 
denote  the  value  of  the  hardening  parameter  at  t^  which 
is  the  same  as  the  value  at  the  end  of  the  last  plastic  pha- 
se preceding  time  t^  and  let  F(t^)  denote  the  loading  func- 
tion of  Eq.  (37).  Then 

< 1 elastic 

F(t. ) if  t.  corresponds  to  an  state  (54) 

=1  1 plastic 
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At  a time  t^+At  , when  the  incident  pressure  is 
Pi ( ti+At ) , w(t^+At)  is  obtained  from  Eq.  (19)  and  then 

l( tj+At)  from  Eqs.  (17)  and  (18).  The  elastic  stress  in- 

crement corresponding  to  the  strain  increment 


Ae  = e(ti  + At)  - e(tL) 


(55) 


is  obtained  from 


( As  ^ = e Ae 


(56) 


in  which  the  subscript  1 denotes  a first  approximation  to 
the  value  which  would  be  obtained  from  Eq.  (28)  if  eP  were 
known.  Correspondinly,  a stress 

s1(ti  + At)  = s(t)  + ( As )x  (57) 

can  be  calculated  and  substituted  into  the  loading  function 


P1(ti  + At) 


< 1 
= 1 


(58a) 

(58b) 

(58c) 


Case  can  occur  only  if  an  elastic  change  occurs.  Thus 


a1(ti  + At)  = s(t  +At) 
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and  the  value  of  all  other  response  functions  at  t=t^+  At 
are  readily  determined. 

Case  b occurs  if  an  elastic  loading  brings  the  stress 

point  to,  or  keeps  it  on,  the  loading  surface  and  Eq.  (59) 

again  holds.  For  both  cases  a and  b,  no  updating  of  the  har- 

* 

dening  parameter  M occurs. 

Case  c indicates  that  plastic  strain  has  occured  and 
that  ( As^  computed  from  (56)  is  an  upper  bound  to  the 
correct  As  . As  a first  correction,  determine  ^ using 
Eq.  (49)  with  and  the  last  update  of  M#  . This 

requires  utilization  of  Eqs.  (43),  (44)  and  (51).  Now 
( A&  is  obtained  from  Eq.  (48)  and  a second  approxima- 
tion to  s(t^  +At)  is  obtained  from  Eq.  (28)  as 

s2(ti  + At)  = E [e  - (aP)1]  (60) 

Using  s2  , a new  value  F^t^  + At)  is  obtained  from  Eq. 
(37).  If  F2  is  ~1, 

F2(ti  + At)  = F(tj  + At)  (61) 

and 


Lu 


a.2(ti  + At)  = + At) 


(62) 


If  ?2  is  still  >1  by  too  much,  a second  iteration  is 
necessary. 


When  Eq.  (6l)  is  held  to  be  valid,  a check  is  made  to 
see  that  plastic  loading  has  occured,  using  Eq.  (45).  if 
Eq.  (45)  is  satisfied  M*  is  updated  using  Eq.  (42). 

The  procedure  outlined  above  is  illustrated  in  Fig.  7. 

As  can  be  seen  in  Fig.  4,  there  are  corners  in  the 
loading  function  along  the  M and  N axes.  At  these  cor- 
ners where  two  loading  surfaces  Fft  and  F^  meet,  as  in 
Fig.  8,  the  strain  increment  is  obtained  from  the  linear 
combination  (Ref.  [13]  ) 


3s 


X 21* 

^ 3s 


(63) 


which  replaces  Eq.  (48).  In  Eq.  (63),  and  are 

obtained  from  Eq.  (49)  with  F = F&  and  F^  , respectively. 
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CHAPTER  4 
NUMERICAL  RESULTS 

A steel  shell  and  water,  with  the  following  geometrical 
and  material  parameters  were  usedi 


1.44  in. 

( . 0366  m) 

20  ft  - 2 in. 

(6.143  a) 

5 ft  - 6 in. 

(1.676 

m) 

42  ft  - 10.5  in. 

(13.07 

m) 

24  ft  - 0 in. 

(7.315 

m) 

15  ft  - 2 in. 

(4.623 

m) 

30  ft  - 11  in. 

(9.423 

m) 

1.94  lb  sec2/ ft** 

(2.073 

x 107  kg/m3) 

14.9  lb  sec2/ ft** 

(15.92 

x 107  kg/m3) 

30  x 106  lb/in2 

(2.068 

x 1011  N/m2) 

4790  ft /sec 

(1460. 

m/s ) 

3600  lb/in2 

(2.48  x 108  N/m2) 

The  above  values  are  the  same  as  for  the  problem  treat- 
ed in  Ref.  [l]  , except  that  the  shell  considered  here  is 
unstiffened. 

Before  proceding  to  the  numerical  solution  of 
this  problem,  the  elastic  portion  of  the  program  to  be 
used,  modified  to  account  for  the  stiffeners,  was  used  to 
solve  the  problem  of  Ref.  [l]  . A very  close  numerical 


check  was  obtained. 
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The  following  numerical  results  were  then  obtained  for 
the  elastic-plastic  problem,  using  a maximum  incident  pres- 
sure of 

(Pl^ax  = 50  psi  (3.45xl05  N/m2) 
and  a time  step 

t ■ .0002  sec  i 

In  Pig.  9»  the  history  of  the  moment  at  the  base  of 
the  shell,  where  yielding  first  occurs,  is  shown  and  com- 
pared with  what  would  have  been  obtained  elastically.  It 
is  seen  that  yielding  first  occurs  at  t **.12  sec,  i.e.  dur- 
ing the  second  spike  of  the  incident  pressure  of  Pig.  3 
(scaled  up  to  a maximum  of  50  psi)  and  that  the  limit  mo- 
ment is  reached  at  t = .22  sec,  i.e.,  during  the  third 

spike  in  Pj  and  at  t a .36  sec  during  the  fifth  spike  . 

It  should  be  noted  that  this  graph  and  those  that  follow 
were  plotted  using  values  for  each  40  time  steps,  i.e.,  for 
increments  of  .008  seconds. 

In  Pig.  10,  the  moment  diagram  when  the  limit  moment 
is  reached  at  the  base  is  displayed,  showing  clearly  that 
the  plastic  zone  due  to  bending  is  confined  to  the  immedia- 
te vicinity  of  the  base. 

In  Fig.  11,  a plot  of  the  value  of  the  loading  function 


at  the  instant  when  the  limit  moment  is  first  reached  is 
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shown.  Note  that  two  internal  plastic  regions  are  primari- 
ly to  hoop  stresses. 

Pig.  12  compares  the  history  of  shell  displacement 
near  the  shell  center  (grid  point  8)  with  what  would  ha- 
ve been  obtained  in  a purely  elastic  response.  After  ini- 
tial yielding  at  t=.128  sec  the  plastic  displacement  ap- 
proaches a damped  oscillation  about  a permanently  deformed 
configuration.  The  displacement  diagram  which  looks  like 
Pig.  13  at  first  yielding,  displays  this  permanent  set  in 
the  displacement  at  grid  point  # = 8 in  Figs.  l4a  and  14b, 
when  its  value  reaches  a maximum  and  a minimum,  respective- 
ly, but  does  not  pass  through  zero. 

Pig.  15  is  a plot  of  the  hoop  force  at  grid  point 
ot  = 8 as  a function  of  time.  The  limiting  value  of  N=NQ, 
is  reached  twice,  at  t = . 24  sec,  ant  at  t = .232  sec. 

At  the  later  time  the  distributions  of  hoop  stress  and  plas 
tic  zones  are  shown  in  Pigs.  16  and  17. 

Pig.  18  illustrates  the  hoop  stress-strain  diagram  for 
grid  point  oc  « 8 . Within  the  global  loading  path  denoted 


by  (p  , some  local  unloading  occurs.  The  corresponding 
moment -curvature  path  for  the  base  of  the  shell  is  shown  in 
Fig.  19. 

Using  the  finite  difference  form  of  the  equation 


the  shear  can  be  obtained.  Its  time  history  at  the  base 
and  distribution  when  the  limiting  moment  is  reached  the- 
re are  displayed  in  Pigs.  20  and  21,  respectively. 

Fig.  22  shows  the  induced  fluid  pressure  on  the  shell 
wall  at  grid  point  «=  8 superimposed  on  the  incident  pres- 
sure there.  The  total  pressure  on  the  wall  is  obtained  by 
superposition  of  these  curves. 
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APPENDIX  I 

In  order  to  avoid  the  possibility  of  violating  Drucker's 
second  postulate,  Eq.  (40),  it  is  not  necessary  to  change 
I*  from  the  expression,  Eq.  (39),  proposed  in  Ref.  [ 2 J , 
to  one  like  Eq.  (36).  Indeed,  in  a problem  involving  a 
general  state  of  stress  in  the  shell,  it  may  not  be  easy  to 
find  a substitute  expression.  Instead,  Eq.  (39)  can  be 
used  and  the  positive  definiteness  of  D (or  the  positive- 
ness of  \ ) checked  by  monitoring  that 


3f 

3s 


(Al) 


If  p>0,  set  =0  in  Eqs.  (49)  and  (53)  and  treat  the 
shell  as  ideally  plastic  for  that  time  increment.  This 
preserves  the  physically  desirable  softening  in  N for 
stress  points  close  to  the  N axis. 

The  calculations  described  in  Section  4 were  repeated 
using  the  modified  procedure  discussed  in  this  Appendix. 
Comparison  of  peak  values  using  the  two  proposed  methods 
are  made  in  Table  I.  It  is  readily  seen  that  both  methods 
produce  essentially  the  same  results. 


A finite  difference  technique  has  been  developed  for 
obtaining  the  axisymmetric  dynamic  elastic-plastic  response 
of  a fluid-filled,  unstiffened  shell  using  shell  theory 
(as  opposed  to  three  dimensional  equations  and  integrating 
through  the  thickness)  in  both  elastic  and  plastic  ranges. 
The  method  is  then  used,  in  two  variations,  to  obtain  the 
response  of  the  wall  of  a water  filled,  nuclear  containment 
vessel  when  the  relief  valve  discharge  piping  is  cleared  and 
fluid  pressures  exerted  on  it  large  enough  to  exceed  the 
elastic  range. 
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TABLE  I 

COMPARISON  OP  RESULTS  OBTAINED  USING 
MODIFIED  METHOD  OF  APPENDIX  AND 
ORIGINAL  METHOD  PROPOSED  IN  SECT.  2.2.2 

* 

ORIGINAL  MODIFIED  DIFF. 


DISPLACEMENTS 

Inches 

Meters 

Inches 

Meters 

(PEAK) 

2.76 

.07 

2.75 

.069 

.40 

1.9 

.0482 

1.89 

.048 

1.26 

3.8 

.0965 

3.79 

.0962 

.20 

3.08 

.0782 

3.02 

.0767 

1.17 

MOMENTS 

Lb . inch 

Joules 

Lb.  inch 

Joules 

(PEAK) 

15493.8 

1750.47 

15503.7 

1751.59 

.06 

-18056.0 

-2039.95 

-17937.4 

-2026.55 

.77 

10936.3 

1235.59 

10895.4 

1230.95 

.40 

- 9363.2 

-1057.84 

- 9359.6 

-1057.44 

. 04 

18195.5 

2055.71 

18010.0 

2034.75 

1.00 

-17094.7 

-1931.34 

-16994.0 

-1919.97 

.60 

13028.4 

1471.94 

12983.2 

1466.83 

.40 

C IMCUMFERENTIAL 


FORCES 

Lbs/inch 

N/m 

Lbs/inch 

N/m 

(PEAK) 

51790.9 

9.1xl06 

51768.3 

9. lxlO6 

.04 

-51743.7 

-9. lxlO6 

-51762.7 

-9. lxlO6 

.04 

37041.9 

6. 5xl06 

36820.2 

6. 4xl06 

.60 

-45762.0 

-8. OxlO6 

-45688.8 

-8. OxlO6 

.20 

51281.1 

8.9xl06 

5H92.1 

8.9xl06 

.17 

-51747.7 

-9. lxlO6 

-51432.5 

-9. OxlO6 

.60 
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Fig.  5.  Course  distribution  of  finite  difference 
molecules  used  for  fluid 
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History  of  shell  displacement  near  center  (grid  pt 


Fig.  13.  Displacement  of  shell  just  after  initial 
yielding  of  base  (t  = .132  sec) 


Fig.  14a.  Shell  displacement  at  t=  .24  sec  Fig.  14b.  Shell  displacement  at  t=  . 31  sec 
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Fig.  15.  Hoop  force  near  center  of  shell 
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Fig.  18.  Hoop  stress  - hoop  strain  diagram  at  grid  point o 


Fig.  19.  Moment  curvature  relation  at  base  of  shell 


Fig.  20.  Shear  history  at  base  of  shell 


History  of  induced  fluid  pressure  (superimposed  on  incident 
fluid  pressure)  on  shell  wall  at  grid  pt.  a = 8 
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APPENDIX  III 


NOTATION 


The  following  symbols  are  used  in  this  paper: 

A = quantity  defined  by  Eq.  (51) 

D = elastic-plastic  tangent  stiffness,  matrix  defined  by 
Eq.  (53) 

E = superscript  denoting  elastic  component. 

= Young's  modulus. 

E = elastic  modulus  matrix,  defined  by  Eq.  (27) 

F = loading  function,  defined  by  Eq.  (37) 

Fa,  = loading  surfaces  meeting  at  a corner.  See  Fig.  8. 
F^ , F^  = successive  iterative  approximations  to  F of  Sec.  3 

Fq  = yield  function,  defined  by  Eq.  (30). 

F^  = limit  function,  defined  by  Eq.  (36). 

G = gradient  of  loading  function,  defined  by  3q.  (43). 

= contribution  to  G of  bending,  defined  by  Eq.  (44) 

IM  = quantity  defined  by  Eq.  (32). 

= quantity  defined  either  by  Eq.  (38)  or  (39). 

1^  = quantity  defined  by  Eq.  (31). 

IMN  = <5uantity  defined  by  Eq . (33). 

L = total  depth  of  model  of  containment  vessel.  See  Fig.  1 


M * subscript  denoting  contribution  of  bending  moment. 

“ longitudinal  bending  moment  per  unit  of  length. 

M*  • hardening  parameter. 

Mq  * yield  moment  in  pure  bending#  given  by  Eq.  (35). 

■ limit  moment  in  pure  bending. 

N = subscript  denoting  contribution  of  hoop  force. 

* hoop  force  per  unit  length. 

N * axial  force  which  causes  yield  in  absence  of  bending, 
given  by  Eq.  (34). 

a,b  = subscripts  denoting  two  loading  surfaces  at  a corner. 
See  Fig.  8. 

c =»  acoustic  velocity  in  water. 

e » total  strain  matrix,  defined  by  Eq.  (22). 

E P 

e , e = elastic  and  plastic  components  of  e,  respectively, 
h * shell  wall  thickness. 

i = generic  subscript  denoting  discrete  time  step. 

k = longitudinal  curvature  of  shell. 

k = value  of  k at  yield  in  pure  bending, 
o 

kE,  kp  = elastic  and  plastic  components  of  k,  respectively. 


I 


depth  of  elastic  shell  wall:  See  Fig.  1. 


i 


5; 


i 


I 


, 


p = dynamic  pressure  in  contained  water. 


Pj  = incident  fluid  pressure  on  shell  wall  caused  by  cleaning 
of  relief  value. 


r = radial  cylindrical  coordinate.  See  Fig.  1. 


ri , r 2 r,  = radii  of  pedestal,  containment  vessel  bottom, 

' and  shell  wall,  respectively.  See  Fig.  1. 

s = stress  matrix,  defined  by  Eq.  (26) . 

s-jy  £2etc.  = stress  matrix  iterations.  See  Sect.  3 

t = time . 

t^  = time  at  i th  discrete  interval, 
u = radial  component  of  fluid  velocity, 
u = (vector)  fluid  velocity, 
v = axial  component  of  fluid  velocity. 

w = (radial)  displacement  of  shell  wall,  assumed  positive 
outward. 

z = axial  cylindrical  coordinate.  See  Fig.  1. 

A = denotes  temporal  increment. 

7 = (vector)  gradient. 

a = angular  discontinuity  in  shell  bottom.  See  Fig.  1. 

6 = quantity  defined  by  Eq.  (Al) . 

5 = denotes  a restricted  temporal  rate.  See  Eqs.  (45) 

and  (47) . 


A 


1 


total  circumferential  strain  in  shell  well. 

P 

fc  = elastic  and  plastic  components  of  t , respectively 
quantity  defined  by  Eq.  (49)  . 

k = value  of  X for  two  loading  function  at  a corner. 
Poisson's  ratio, 
uniaxial  yield  stress, 
density  of  fluid. 


density  of  shell  wall. 
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